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1 Introduction 



Neutrinos with magnetic or transition magnetic moments (//) propagating in transverse mag- 
netic fields (Bj_) will undergo spin precession |l], ||, |3j or spin-flavour precession |4], ||, |5|, ||, 
respectively. As is well known, both types of precession can have important astrophysical 
implications. In the presence of matter, the neutrino spin-flavour precession can be reso- 
nantly enhanced [[5], ||. The resonant spin-flavour precession (RSFP) is analogous to the 
MSW effect in neutrino oscillations |7|, ||. If neutrino mixing exists, an interplay between 
the RSFP and the MSW effect is possible. In particular, the matter-enhanced neutrino 
oscillations can assist the RSFP in giving rise to time variations of the solar neutrino flux 
by improving the adiabaticity of the latter [[J. 

In most of the previous studies of the spin and spin-flavour precessions of neutrinos in 
transverse magnetic fields it was assumed that the direction of the field is fixed. However, 
in the physically interesting cases of three-dimensional field configurations (e.g., in the Sun, 
supernovae, etc.), neutrinos traversing the magnetic field will feel, in general, a field whose 



direction is changing along the neutrino trajectory [T^, 11, 12, 13, 14]. The presence of such 
twisting magnetic fields can lead to new interesting phenomena. In particular, the neutrino 
spin precession can be resonantly enhanced, the energy dependence of the probability of the 
resonant spin-flavour precession can be significantly changed, the lower limit on the value 
of fiB± required to account for the solar neutrino problem can be slightly relaxed and, in 



addition, new realizations of RSFP become possible |12. 13|. 

The RSFP of neutrinos in twisting magnetic fields was considered so far under the sim- 
plifying assumption of absence of neutrino mixing in vacuum. In the present paper we study 
the general case: we assume the existence not only of flavour-off-diagonal neutrino magnetic 
moments, but also of flavour neutrino mixing in vacuum. This is a natural possibility if no 
additive lepton number is conserved. Under the above condition both resonant oscillations 
and RSFP of neutrinos are allowed. We show that if they take place in twisting magnetic 
fields, qualitatively new phenomena become possible: resonant enhancement of the helicity- 
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flipping ALi = 2 neutrino-antineutrino conversions un <-> vir (I = e,//, r), and merging 
of three different transition resonances. Moreover, for a definite sign of the angular veloc- 
ity of the magnetic field rotation, the order of the RSFP and MSW resonances in matter 
with monotonically decreasing density can change. All these effects can have an important 
impact on solar and supernova neutrinos. 



2 Neutrino Evolution Equations 

Consider a system of two flavour neutrinos and their antiparticles, say u eL , P eR , v^ L , and u^r. 
We shall assume that v e - mixing is generated by a Majorana neutrino mass term, and 
that there exists a v e L - v^r transition magnetic moment /i. The equations for the flavour 
neutrino wave functions, describing the propagation of neutrinos in matter and transverse 
magnetic field in the case of interest, can be written in the form: 

s 2 5 vB ± (t)e i(t> W \ I u eL y 

-nB^e-W) S2 5 | P eR 

r(t)N(t) + c 2 5 | v^ L 

r(t)N(t) +c 2 5 ) { v^r , 

(1) 

Here the angle <j)(t) defines the direction of the magnetic field Bj_(t) in the plane orthogonal 
to the neutrino momentum, B±(t) = |Bj_(i)|, 

N = V2G F (n e -n n /2), r = — ^ , 5 = Am 2 /4E, 

s 2 = sin2# , c 2 = cos 26*o, (2) 

where Gf is the Fermi constant, n e and n n are the electron and neutron number densities, 
E is the neutrino energy, Am 2 = m\ — m\, mi and m 2 being the masses of two Majorana 
neutrinos and # being the neutrino mixing angle in vacuum. The zeros in the effective 
Hamiltonian in eq. (1) are related to the fact that diagonal magnetic moments of Majorana 
neutrinos are precluded by CPT invariance. 
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It is instructive to consider the evolution of the neutrino system in the reference frame 
which rotates with the same angular velocity as the transverse magnetic field The 



corresponding transformation of the neutrino wave functions reads: v = Sv' , where v is the 
column vector of state of the neutrino system in eq. (1), and S = diag(e l 2 ; e~*2 , e 1 ? , e~ l 2 ). 
Since S is a diagonal matrix, the above transformation does not change the probabilities 
of the neutrino transitions and therefore we will keep the same notations for the neutrino 
wave functions in the new frame. In the rotating frame the phase factors disappear from 
the nondiagonal terms (see eq. (1)) and the evolution equations take the form: 
/ 
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(3) 



where = d<p/dt. 

The 4x4 matrix on the right-hand side of the Schroedinger-like equation (3) defines the 
time-dependent effective Hamiltonian of the system, H. Note that the effect of rotation of 
the magnetic field amounts to the appearance of additional terms of the same absolute value 
in the diagonal elements of H . These terms have the same signs for the neutrinos possessing 
the same helicity and opposite signs for the neutrinos of opposite helicities. The indicated 



features have a simple physical interpretation ]nj. According to classical mechanics, a 
system with angular momentum s acquires an additional energy AE = — (sf2) upon a 
transformation to a reference frame rotating with an angular velocity O. In our case Q = <fi, 
s = 1/2, and the projections of s on ft are equal to —1/2 for the neutrinos, and to +1/2 for 
the antineutrinos. 

Let us consider the resonant properties of the neutrino system under study. The diagonal 
elements of the Hamiltonian in eq. (3) define the energies of the "flavour levels", i.e. of 
VeL-, ^eR, ^l, and v^r, in the absence of vacuum mixing and magnetic field; we shall denote 
them respectively as H e , H s , H^, and H^. Resonances take place at densities at which 



the splittings between the eigenvalues of the effective Hamiltonian H with s 2 5 ^ and 
lxB± ^ are minimal. The resonance densities can be approximately found by equating 
pairs of diagonal elements of H: H a (n e ,n n ) = Hp(n e ,n n ), a ^ (3, a, (3 = e, e, //, Jx (i.e. they 
correspond approximately to the values of N at which the flavour levels cross). 

With four neutrino states there are altogether six resonance conditions, and in what 
follows we will consider the nature and the properties of the corresponding resonances. 

Two conditions, H e = and H E = H^, can be written as: 

/\ ?T?,^ 

V2G F n e = — — cos 26» [y eL <-» v^), (4) 
/\ th 2 

V2G F n e = — — - cos26» (i> eR <-> v^ R ). (5) 

They concern the transitions (given in parentheses) between neutrinos having the same 
helicity but different flavour, and thus correspond to the MSW resonances. The condition 
for the neutrinos (antineutrinos) can be satisfied for Am 2 cos26> > (Am 2 cos2# < 0). 
Eqs. (4) and (5) do not depend on since the flavour neutrino states involved have the 
same helicity and the magnetic field rotation influences both states equally. Magnetic fields 
(twisting or nontwisting) do not modify the MSW resonances. 

The next two conditions, H e = and H s = H^, have the form: 

/\ 777 , 

V2G F (n e - n n ) = cos26» - [v eL <-»■ v^r), (6) 

V2G F (n e -n n ) = — — cos26> o -0 (v eR ^v^ L ). (7) 

The neutrinos involved in the relevant transitions have different flavours, as well as opposite 
helicities. Thus, eqs. (6) and (7) define the resonance conditions for the spin-flavour preces- 
sion: eq. (6) refers to the channel v eL <-> u^ R , while eq. (7) refers to the CP-conjugated one. 
Since the neutrinos involved have different helicities, the magnetic field rotation modifies 
the resonance conditions: both eq. (6) and eq. (7) depend on 0. In a medium with (n e — n n ) 
of a definite sign (like the matter of the Sun) and = only one of these two conditions 
can be satisfied depending on the sign of Am 2 cos2# - In contrast, in the case of twisting 



magnetic field both conditions can be obeyed simultaneously. This can happen for < 
(0 > 0) when n e > n n (n e < n n ) provided |0| > \(Am 2 /2E) cos26 , |. Obviously, the two 



conditions will be fulfilled at different densities. 



Finally, equating the elements H e = H E and i? M = one finds: 



V2G F (n e -n n /2) 
V2G F (-n n /2) 



0/2 (z/ M <-»■ 



0/2 (z/ e <-> z/ e ), 



(8) 



(9) 



Equations (8) and (9) represent the resonant conditions for the channels u eL <-> z/ efi and 
*Vr, respectively. In these transitions the helicity is flipped and the electron or the 
muon lepton number is changed by two units. Thus, we deal here with a qualitatively new 
phenomenon: neutrino-antineutrino helicity flipping ALi = 2 resonance conversions. 

In a medium with nonvanishing (n e — n n /2) and n n (like the Sun's matter) conditions 
(8) and (9) cannot hold in a nonrotating magnetic field. For fixed direction of the field 
rotation only one of these two conditions can be fulfilled: if, for example, n e > n n /2, for 
< the resonance takes place in the electron neutrino channel, v eL <-> P eR , whereas for 
> it occurs in the muon neutrino channel, u^ L <-> u^ R . 

It should be emphasized that the neutrino-antineutrino transitions are forbidden in the 
lowest order of the interactions described by the effective Hamiltonian in eq. (3): the matrix 
elements of H which mix neutrino and antineutrino states of the same flavour are zero. As 
we have mentioned earlier, this follows from the fact that Majorana neutrinos cannot have 
nonzero diagonal magnetic moment. In Section 4 we will show that the mixing between 
neutrinos and antineutrinos is induced in the second order of the perturbation theory in the 
interactions described by H . 
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3 Merging and Permutation of Resonances 



The energies of the flavour neutrino levels, H a , have different <fi dependence. Consequently, 
the resonance densities associated with the various neutrino transitions are different func- 
tions of 0. Therefore the magnetic field rotation modifies the level crossing scheme. In 
particular, one can expect such effects as reduction (or increase) of the spatial distance 
between the resonances, or merging of resonances, or permutation of resonances, to take 
place. 

In what follows we shall consider the dependence of the level crossing scheme on the 
value of 0. We shall assume for simplicity that the chemical composition of the medium 
does not change appreciably along the neutrino path, so that r(t) = const. In this case the 
level energies H a depend on one density parameter only: N(t). It can take both positive 
and negative values; in particular, N(t) is always positive in the Sun, but can be negative 
in a strongly neutronized medium like the one existing in the central regions of collapsing 
stars. 

It proves convenient to subtract the matrix H e -I, proportional to the unit matrix /, from 
the Hamiltonian H in eq. (3). This is equivalent to a multiplication of all wave functions by 
the same phase factor exp(—iH e t) which, obviously, does not change the neutrino transition 
probabilities. The diagonal elements of the Hamiltonian (if — H e ■ I) read: 

{H e ,H s ,H^H p } = {0, -2N-j>, -(l + r)N + 2c 2 5, -(1 - r)N + 2c 2 5 - 0}. (10) 

The resonance values of the density parameter N, corresponding to the different neutrino 
transitions, can be easily deduced from eq. (10). For flavour neutrino transitions one has: 

A7<;> = (") 

where the plus sign corresponds to the u eL <-> u^ L channel, and the minus sign, to the CP- 
conjugated one v e R <-> v^ R . The resonance values of N for the spin-flavour conversions are 



6 



equal to 



±2c?5 - 



N sf(sl) = 1 _ r \ (12) 



where the plus (minus) sign corresponds to the v eL <-> u^ R (u eR <->• v^ lL ) transitions. Finally, 
the resonances in the neutrino-antineutrino transitions v e L —>■ v eR and u^l — > take place 
respectively for 

^ee-=-|, (13) 

^ - (14) 

We shall assume for simplicity that <fi = const along the neutrino trajectory. Then 
the energy levels H a = H a (N) are straight lines (see Figs. 1 and 2). After the phase 
transformation of the neutrino wave functions i) H e = by definition, and ii) the muon 
neutrino energy level H^ does not depend on 0. The antineutrino energy levels H s and 
Hp have the same <fi dependence. Consequently, when varies, H e and do not change, 
whereas H g and Hp are shifted by the same amount. Finally, note that for r < 1 (which is 
realized in the Sun) the strongest dependence on N is exhibited by H s , the dependence of 
Hp being somewhat weaker, while that of Hp is the weakest. 

We shall assume that neutrinos are propagating in matter with density decreasing mono- 
tonically along the neutrino trajectory. Consider first the possibility C2<5 > 0. Let us recall 
the level crossing scheme in the case of nonrotating magnetic field. With = (Fig. 1) 
there are only two resonances for N > 0: the spin-flavour, z/ e ^ <-» v^r, and the flavour (i.e. 
MSW) one, v e ^ «-> u^l- The spin-flavour resonance is located at higher density as compared 
to the density of the MSW resonance || |TB[ : 



N sf _ 1 + r 



> 1. (15) 



N f 1 - r 

Therefore if v e L is produced at a density Nq > N s f and the adiabaticity condition is fulfilled 
at the spin-flavour resonance, v e ^ will be converted into (provided the resonances do 
not overlap). If Nq < N s f, or if the adiabaticity at the spin-flavour resonance is strongly 



violated, v e L will be transformed completely or partially into v^l, depending on the degree 
of adiabaticity at the flavour resonance (N = Nf). 

Suppose now that 7^ 0. The effects caused by being nonzero depend on the direction 
of the field rotation and we shall discuss first the case < 0. With |0| increasing, the u^r 
and u eR energy levels become higher and the following changes take place. 

a) < |0| < 2c2# (Fig. 2a). A u e L — u eR level crossing appears for iV > 0, so that under 
certain conditions (see Section 4) u eL can undergo a resonant conversion into u eR . 

b) 2c 2 5 < |0| < 4c 2 5/(l +r) (Fig. 2b). For N > a resonance appears in the — v eR 
channel. This means that two different types of spin-flavour conversion can take place in 
a medium with N > 0: u e L <-> P^ R , and u^l ^ v eR . Moreover, in such a medium u e L 
can be transformed into v eR via two consecutive resonant conversions: first the MSW one 
v e L — > VjiL (at iV = Nf) and then the spin-flavour one — * v eR (at = N s y). We shall 
call this process a two-step u eL —> u eR resonant conversion. 

c) At = 0i m , where 

X 4026 (1R\ 

01m = (16) 

(Fig. 2c), three resonances merge in one point: Nf = N s j = N eS = 2c2<5/(l + r ) ( see ec L s - 
(11)-(13)). In particular, the two resonances present in the two-step u eL — > u eR transition 
merge and a strong v eL — > v eR conversion can take place in one resonant region (see Section 
5). 

d) For |0| > 1 0i m I (Fig. 2d) the flavour and the spin-flavour resonances permute and 
if iVo < N s f, a v eL propagating from high to low densities will cross first the u eL <-> u eR 
resonance and later the v eL <-> u^ L one. If the adiabaticity condition is fulfilled, u eL will 
undergo efficient conversion into u eR in the direct v e ^ — > v eR transition before reaching the 
MSW resonance (for which the adiabaticity parameter can be much bigger). 

For large values of |0| the energy level Hp is far away from the energy domain where the 
v e L, v iiL and u^ R levels cross, so that in the corresponding density region the v^ R decouples 
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from the dynamics of the other three neutrinos (see further, Section 5). 

For positive values of the antineutrino energy levels H s and become lower as 
increases. The following specific features of the level crossing scheme are worth noticing. 

e) If 4> > 0, a resonance in the <-> v^r channel appears for iV > (Fig. 2e). If the 
adiabaticity condition is fulfilled for this resonance, v e L can undergo a transition into 
through the two-step resonant process of v e L — > v^r and v^r — > v^l, imitating the flavour 
conversion v eL — > u^ L . 

f) For = 2m , where 

2m = 4c 2 6-^— , (17) 
1 + r 

(Fig. 2f) the resonances in the three channels u & l <-> z/ m l, z/^r <-> tyj,, and z/ M ^ <-> z/ e £ merge 
in the point N — Nf. 

g) When > 2rrt , permutation of the flavour and the spin-flavour resonances occurs, 
so that now Nf > N s f in contrast with the case of nontwisting magnetic field. Therefore, if 
a u e L moves from sufficiently high to low density regions and the adiabaticity condition is 
fulfilled in the flavour resonance, it will be transformed into u^l rather than v^r. 

h) For > 2c 2 5 (Fig. 2h) the spin-flavour resonance disappears if N > (in this case 
N s f < 0). Thus, for N > only two resonances exist: in the v e L <-> u^l and u^l <-> v^r 
transitions, the latter occurring at larger density. For sufficiently large positive values of 
the v e R state decouples from the rest of the system. 

Let us consider next briefly the transitions of V eR . The most interesting effect is the 
v e R — * ^eL conversion which can take place for negative values of (see Figs. 2a-2d). It can 
proceed as a direct resonant conversion (Figs. 2a,b), or as a two-step process v e R — > v^l 
v e L (Fig. 2d). In addition, v e R —> v e L conversion is possible in the case of merging of three 
resonances (Fig. 2c). This may have important implications for the neutrinos emitted in 
collapsing stars. 

Finally, let us discuss the interesting special case of r = 1. As it follows from eq. (2), it 
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corresponds to propagation of neutrinos in an isotopically neutral medium, n e = n p = n n , 
and can be realized in wide density regions in the supernovae. For r = 1 the matter density 
term (N) no longer enters the resonant conditions for the spin-flavour conversions (eqs. (6) 
and (7)), which take the form: 

= ±2c 2 5. (18) 

Thus, the spin-flavour conversions in the case under discussion can take place only in rotating 
magnetic fields. Further, for r = 1 the dependences on iV of the u e R and u^l energy levels, 
as well as those of the u eL and u^ R energy levels, are the same: the energy levels H a as 
functions of N form two pairs of parallel lines (Fig. 3). For = —<j) m = —2c 2 5, merging 
of two pairs of resonances, namely, of those in the channels v e ^ <-> and v e ^ <-> P c r, and 
in the channels P^r <-> and P e R <-> P^r, takes place. Moreover, the levels of u^l and 
v e R coincide now. This means that the v e R <-> precession will proceed with practically 
maximal amplitude for any value of N, except in the vicinity of the triple resonance points. 
For = +2c 2 5 the levels of P^r and u e L merge, and now the u e L <-> P^r precession will 
proceed with maximal amplitude. 

One can perform similar analysis in the case c 2 5 < 0. The results can be obtained 
formally from those derived above for c 2 5 > (see a) - h), etc.) by replacing c 2 5 by 
—c 2 5, and by interchanging uil and P\r. For instance, if < there will be again a 
resonance in the transition u eL <-> P eR , and for = lrn = 4c 2 S/(l + r) three resonances, 
this time in the channels v eR <-> P^r, u eL <-> P eR and u eL <-> P^r, will merge in one point 
(N = N f = N eS = N sf = -2c2(J/(l+r)). 

As we see, an electron neutrino propagating in matter and magnetic field can be trans- 
formed into any other state of the neutrino system. The neutrino final state will represent, in 
general, a mixture of the states of the four neutrinos: is e L, P € r, v^l, and P^r. Their relative 
admixtures in the final state are determined by i) the density N at which the neutrinos are 
produced, ii) the order of the different resonances crossed by the neutrinos, iii) the degree of 
overlapping of different resonances, and iv) by the degree to which the adiabaticity condition 
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is fulfilled in each of the resonances. Magnetic field rotation can change the order of the 
resonances and also can indirectly influence the relevant adiabaticity conditions. In the case 
of nonuniform field rotation <ft enters the adiabaticity conditions directly [|12|, Magnetic 
field rotation can change significantly the composition of the final neutrino state. The most 
interesting new effects from a phenomenological point of view are the v e L <-> v e R transitions 
and the merging of different resonances; we will discuss these effects quantitatively in the 
next two sections. 



4 v — v Resonances 

Let us consider the v e L <-> v £ r transitions in the vicinity of the v e L — v e R level crossing. As 
we have mentioned earlier, to lowest order of the perturbation theory in the interactions 
described by H there is no v e L — u e R mixing. The latter arises as a result of the interplay 
of the two types of mixing present in H in eq. (3): the flavour mixing and the mixing 
generated by the magnetic moment interaction with the magnetic field. More specifically, 
according to eq. (3) the transition v e L — > 9 e R can proceed in two ways: 

VeL — > v^l — > u eR , (19) 

V e L > ViiR > V eR , (20) 

i.e. via v^l or u^r in the intermediate state. In (19) the first transition is due to the flavour 
mixing (S2S) and the second is induced by the interaction with the magnetic field (fiB±), 
and vice versa in (20). Since the system of evolution equations (3) formally coincides with 
the Schroedinger equation, one can make use of the perturbation theory to get the following 
amplitudes of the processes (19) and (20): 

M 1 = -s 2 5fiB ± P(v flL ), (21) 
M 2 = fiB ± s 2 5P(^ R ). (22) 
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Here and P{i>hr) are the propagators of the system with [v^r] in the intermediate 

state: 

PM^ R )) = (H e -H m )- 1 . (23) 

The opposite signs in eqs. (21) and (22) are related to the fact that the matrix of Majo- 
rana magnetic moments of neutrinos is antisymmetric. Evidently, for Piy^ = -P(z/ m _r) a 
complete compensation between the two matrix elements takes place and the v e L — > v e R 
transition is not possible even in the second order of the perturbation theory. Since the self- 
energies of v^i, and v^r are different in matter, the cancellation between the two amplitudes 
(21) and (22) is not complete. Direct calculation in the second order of perturbation theory 
yields the following amplitude of the v eL — > v eR transition: 

2s 2 5 uBArN - 6/2) , , 

M(u eL - u eR ) My • M 2 2 f J K V- L — 24 

[2c 2 5 - (1 + r)N][2c 2 5 - 6 - (1 - r)N] 

Eq. (24) defines an effective coupling which mixes v e L and i> eR , or in other words, a nondi- 
agonal v e ^ — v e R term in the higher order effective Hamiltonian of the neutrino system. 

From eq. (24) there follow several qualitative conclusions, i) The v e L — v e R mixing 
can exist even in the absence of magnetic field rotation [B, [Tj|. ii) Since rN ~ n n , in a 



nonrotating field the v e ^ — v e R mixing vanishes when the neutron concentration is zero [|15 

iii) For sufficiently small n n , the v e L — v e R mixing is induced by the magnetic field rotation. 

iv) The v e L — v e R mixing is a function of 6, and depending on the sign of 6 (i-e. on the 
direction of the field rotation) a twisting magnetic field can either increase or reduce the 
v e L — v e R mixing, and consequently can enhance or suppress the probability of the v e h — > v eR 
transition. In particular, in the u^l — v^r resonance point, where 6 — 2rN, one has M = 0. 

In order to study the properties of the v e h — > v eR resonant conversion one has to solve, 
in general, the system of neutrino evolution equations (3). However, in some physically 
interesting cases the v e h — v e R system decouples from and v^r and its evolution can be 
analyzed separately. This is realized when the v e L ^ v e R conversion resonance is sufficiently 
well separated both in density and in energy from the resonances associated with the v e L <-» 
v^l and VeL v^r transitions. In other words, the v eL — v e R resonance should not overlap 
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with the v e L — VjiL and v e L — v^r resonances, and the mixing of v e L and v e R with and 
i>nR should be sufficiently small. 

There are altogether four general conditions ensuring the decoupling of the v e i — v e R 
system from and v^r. - H e \ > 2|# Me |, \H n - H e \ > 2|if Ae |, \H^ - H E \ > 2\H^\, 
and \Hp — H s \ S> 2\Hj ie -\. Using the effective Hamiltonian (if — H e ■ I), one obtains in the 
"vicinity" of the v e L — v e R resonance, where \c 2 5\ \N + 4>/2\, 

|^|,|^| > \2s 2 5\,2fiB ± . (25) 

The conditions for nonoverlapping of the resonances in density scale can be written as 

\N e e - N f(f) | > 2 tan 26 N f(f) , (26) 
\N eg - N sf{s - f) \ > 2 (fiB ± /\c 2 S ± <P/2\)N sf(s - f) , (27) 

where tan26*o Nf/f\ and ({iB±/\c 2 5±<fi/2\) N s j^, are essentially the half-widths of the MSW 
and the spin-flavour resonances, the width of the v e L — v e R resonance being considerably 
smaller. It is easy to show that (26) and (27) are always satisfied when inequalities (25) 
hold. In a medium with n e > n n conditions (26) and (27) can be realized for <fi < and, e.g. 
c 2 5 > if, for instance, 2c 2 5 <C \<fi\, or if 2c 2 5 ^> \<fi\. In the first case for r = 1 we have (see 
Fig. 2d) Nf <C N eS "C N 8 f, while in the second the inequalities N eS <C Nf, N s f (see Fig. 2a) 
hold. If in the second case, e.g., Am 2 is sufficiently large, only the resonance condition (8) 
can be fulfilled, and so only the v e L — > v e R conversion can be resonantly enhanced. 

If conditions (25) hold, one can block-diagonalize the evolution matrix in eq. (3). The 
resulting 2x2 submatrix governing the evolution of the v e i — v e R system [] takes the form: 



' 2fiB ± s 2 6(3 ^ 

v 2fiB± s 2 5 /3 -2N - j> + 2(s 2 5 2 - /j 2 B 2 ± )(3 ) 



{21 



1 In fact, this matrix describes the evolution of v' e and v' e states obtained as a result of block- 
diagonalization of the initial Hamiltonian, but they coincide with v e and v e up to small corrections of 
the order of 2^^ 2 B\ + s 2 2 5 2 /\ - rN + 4>/2 ± 2c 2 S\. 
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The nondiagonal term, H eS , in the evolution matrix (28) f\ coincides with the effective 
v e L — v e R mixing parameter obtained earlier in the second order of the perturbation theory, 
eq. (24). In the approximation used the parameter j3 can be written as 

P . ~ rN + * /2 . (29) 

The resonance condition following from (28) is 

N + <j)/2 = (s\5 2 - /i 2 5|)/5. (30) 

It is more accurate than the previously found condition, eq. (8), but the correction is 
small. For fiBj_ = s 2 5 the two conditions coincide. The — u e R mixing angle in matter, 
corresponding to (28) is given by: 

tan20 = ■ 2 -^l - (31) 

AT + 0/2- (s 2 2 5 2 - /J 2 B 2 ± )(3 V ; 

Let us consider now in greater detail the nondiagonal element H eS of the v e ^ — v eR 
evolution matrix (28). This element determines the level splitting at the v e L — ^eR resonance, 
as well as the degree of adiabaticity of the u eL <-> u eR transition. It can be represented as 

H ee = -s 2 6»B ± H »~* Ili , (32) 
tip tip, 

and it follows from eq. (32) that under the conditions ensuring the validity of the evolution 
matrix (28), the v e L — v e R mixing is always smaller than each of the generic first order mixings 
s 2 S and /i-Bj_. If one of the two elements and is much bigger than the other, H eS will 
be suppressed by the smaller of the two: e.g., if \Hp\ ^> we have H eS ~ — s 2 5 fiB±/H^. 
At the u e L <-> v e R resonance point (where A" = Aeg = —0/2) H eS can be written as: 



^ - -{2c 2 5Y-(l+rfN% -l-{N e ,/N m f 

? llR 2(l + r)/(l-r) 2 

S2d ^ N eB (l - N sf /N eg )(l - N.J/N*) ' {66) 



2 We have again subtracted a matrix proportional to the unit matrix in such a way that the upper diagonal 
element in the v e L — v e R evolution matrix (28) thus derived is zero. 
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where ti = tan 26 , and Nffj\ and -^/(s/) are the MSW and the spin-flavour resonance 
densities determined in eqs. (11) and (12). It follows from eq. (33) that the v e i — v e R 
mixing term at the resonance increases when N eS approaches Nf, i.e when the v e L — v e R 
resonance approaches the MSW resonance []. However, in the approximation we use, which 
implies a sufficiently good separation between the different resonances, N eg cannot be very 
close to Nf. Using criteria (25) and choosing N eS = Nf ± 2t2 Nf in eq. (33) one obtains for 
the allowed maximal value of H P 



L eei 



B?*w H -<x±kr (34) 

and for small values of t-i (2 1^1 ^ 1) one finds: \H^ S \ = j/i-Bj_. Thus, the "maximal" value 
of | H^ s | in the case considered is only by a factor of 4 smaller than the spin-flavour mixing 
element in H. 

From (33) it follows that \HH S \ increases also when N eS tends to N s y, i.e. when the 
v eL ~ v e R resonance approaches the v eR — resonance. The "maximal" value of \S^ S \ 
compatible with constraint (27) in this case is ||s2<5| and is reached for yU-Bj_/(c2<5) <C 1. 

The mixing parameter Hl§ s enters the expression for the adiabaticity parameter char- 
acterizing the v e L <r-¥ v eR transitions: = (2ffJ| s ) 2 /Hg^ res , where H s ^es is the derivative of 
H E at the resonance, H s being the v eL — v eR level splitting (see eq. (10)). Using the lowest 
order expression for H s , on finds: 



L ee - 2 (Nee/Nf) 



2 



2 ^ ' N es [l-(N eS /Nf) 2 f 1 ' 

Here L eg = N eS /\(N + <p/2) res \, where the value of (N + <p/2) is taken at the resonance point. 
For ^> 1 the v e L —>■ v e R transition is adiabatic and v e L can be completely converted 
into u e R. Again, when N eE approaches Nf (and/or N s j) the adiabaticity of the v eL — > v eR 
transition improves. 



3 We assume here that C2S > 0. 
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5 Triple Resonance. Neutrino - Antineutrino Oscilla- 
tions. 



We shall analyze here the case when resonances merge, which is opposite to the one discussed 
in Section 4. The mutual influence of the different neutrino transitions becomes very strong 
in this case. 

If = 0i m < (see Section 3, eq. (16)), three level crossing points, namely, the 
v eL — v^l-, v e R — v^l, and u e R — v e L, merge at N = Nf. The condition of merging depends on 
the neutrino energy E; for fixed E the resonances merge at definite value of <fi. If neutrinos 
possess continuous energy spectrum, the merging will occur for a given <fi < at definite 
value of the energy, E m : 

(1 + r)<p 

Thus, the specific phenomena which take place at the merging point will have a resonance 
character as functions of energy. 

Suppose = 0i m and consider the dynamics of the neutrino system in the vicinity of 
the merging point which we shall call also point of the triple resonance. Now three neutrino 
states are participating in the transitions and only one state, that of P^r, can decouple in 
the triple resonance region. One has for this state 

H f = - (1 -r)N + 2 ^±ll (37) 

If S2S,/iB± <^ Hp (which implies that s 2 1 and fj,B± 4c2<5), it is possible to perform a 
3 — 1 block-diagonalization of the neutrino evolution matrix in the triple resonance region. 
This leads to the following effective Hamiltonian for the {v eL , v eR) v^l) system: 

( -Ugt -e s 2 S \ 



H' 



-e -2N + ^-^ff -fiB ± 
y s 2 5 -fiB ± -N(l + r) + 2c 2 5 J 



(38) 
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where 

e = if^litM (39) 
Hp 

corresponds to a small direct u eL — u eR mixing of the second order, generated by the prop- 
agation only of V^ R . As we shall see, the mixing induced by the z/^-exchange is resonantly 
enhanced and turns out to be of the order of the generic mixings. The diagonal elements 
also acquire corrections due to the propagation of u^r which are written explicitly in (38). 
They modify condition (36) and the density at the triple crossing point by terms ~ e. 

Let us find the level splitting and the mixing in the neutrino system described by the 
Hamiltonian (38), in the merging point. The eigenvalues of the Hamiltonian (38) read: 

H^f-^H^-f-^H,^ e'. (40) 

Here 



/ = vW) 2 + M) 2 (4i) 

and 

e' = -esin2cj, (42) 

where 

tanc<j = — — (43) 

is the ratio of the leading order flavour and spin-flavour mixing terms. The parameter / 
determines the level splitting in the leading order (i.e. neglecting the terms ~ e): H 1 —H 2 = 
2f, Hi — H 3 = H 3 — H 2 = f . Moreover, the flavour and the spin-flavour mixing terms enter 
/ symmetrically. (For comparison, notice that the splitting in the flavour and in the spin- 
flavour resonance points are equal, respectively, to 2s 2 5 and 2fiB±.) The H 1 — H 2 splitting 
does not depend on the z/^-level correction e', while H 3 shifts towards Hi or H 2 depending 
on the sign of e'. 

The orthogonal matrix diagonalizing the Hamiltonian (38) is given to leading (first) 
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order correction in e/ f by 



Sm — 



\ V2 



7# + 7 
775 ~~ £ 



sing; 

_ COS_t£ 



7 



COSCJ 

sincj 



-73 -f -7 cos 2a; 



.f 



(44) 



where the parameters a, 7 and £ are all of the order of e//: 

1 e n 3 • 2 x 
a = — = — cos a; ( 1 sm u ) , 



7 



1 6 • ft 3 2 V 

-— == — smw 1 cos a>), 

V2f y 2 ; ' 



1 e 

— sin 2a;. 



(45) 

(46) 
(47) 



S 4^2/ 

Consider now the properties of the various neutrino transitions which can take place 
at the merging point. If the density and the chemical composition of the medium, and 
the magnetic field (i.e. N, r, B± and 0) do not change along the neutrino trajectory, 
the propagation of neutrinos has a character of pure oscillations-precessions with constant 
depths and periods. There are altogether three transition channels, v e ^ <-> v^, P eR <-> v^l, 
v eL ^ v eRl and we shall concentrate on the qualitatively new mode u eL <-> P eR . Using (40) 
and (44) we find for the relevant oscillation-precession probability in a medium with N, B±, 
and having values satisfying the merging condition: 

1 3 

P(v e L — ► Veil) = sin 2 2a; sin 4 - ft + sin 2 2a; sin 2 -e't cos ft + d sin e't sin ft, (48) 

where 

sin 2 2a; ^li^ll (49) 

and d = 0(e). As we have already emphasized, the z/ e L <-> u eR oscillations-precessions can 
exist only in the presence of both vacuum mixing and magnetic moment interaction. The 
first term in (48) corresponds to the lowest order, obtained when the term e in (38) is 
neglected. The depth of the oscillations, sin 2 2a;, does not exhibit any suppression related 
to the higher order direct u eL — u eR mixing. Moreover, in the symmetric case when the 
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vacuum mixing is equal to that induced by the magnetic moment interaction, s 2 5 = fiB±, 
the oscillation depth becomes maximal. It should be stressed that this is possible only in 
twisting magnetic fields. For s 2 5 7^ /jlB± the oscillation depth is less than unity and it 
decreases when the difference between s 2 5 and fiB± increases. 

The v eL <-> u eR oscillation-precession length is given by 

l. s * 2 " (50) 

f \/M) 2 + (pB^ 



For (s 2 S) <C (fxB±), this length l eg is two times bigger than the usual precession length, l p : 
l eg — 2/ p = 2n/(fxB±). For (s 2 S) 3> (/i-B_i_) it is twice as big as the oscillation length in 
matter at resonant density, l osc : l ee - = 2l osc = 87iE/(Am 2 sin2# ). In the symmetric case 
when (s 2 S) = (fJ,B±) one has: l ee - = \[2l osc = \f2l p . Finally, note that P(u e L — > v eFt ) depends 
on the fourth power of sin(7rt// e g) rather than on the second power, as in the usual case. 

These features of the v e L v e R oscillations are consequences of the facts that there are 
three levels involved and that the splitting between the levels is determined by /. More 
precisely, according to (44), v e h and u e R are (to leading order) orthogonal combinations of 
two states: uq with energy H 3 = 0, and v' which itself is a maximal mixture of two states 
with opposite energies ±f. The amplitudes of the v — * z/ and v' — > v' transitions are 
equal, respectively, to 1 and cos ft. Their sum in the amplitude of the v e L — * v £ r transition 
gives (1 — cos ft) = 2 sin 2 ^ft. As a consequence, the v e L — > v e R oscillation probability is 
proportional to sin 4 \ ft. 

The second term in (48) is generated by the direct v eL — u e R mixing in (38). It gives 
a long-period modulation of the oscillation probability. For the corresponding modulation 
length, l mo d, one finds: l mo d = 47r/ (3e') 3> l e g- The first two terms in (48) can also be written 
as I sin 2 2u (1 + cos 2 ft — 2 cos ft cos fe't), which implies that the modulation leads to the 
oscillation depth varying between sin 2 2u and \ sin 2 2uj. 

Let us turn next to the neutrino transitions in the region of the merging point in the case 
of varying matter density. The efficiencies of the resonant conversions depend on the degree 
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of their adiabaticity. In what follows we will neglect the corrections ~ e due to the v^r level. 
This can be done, in particular, because the time of propagation in the resonance region 
(t ~ Ar res ) is much smaller than the period ~ 1/e. In addition, as we have established 
above, the relevant level splittings are determined by s 2 S and /iB±, the terms ~ e giving 
only small corrections. 

Consider first the dependence of the energy levels on N (Fig. 4). Its specific features 
are determined essentially by the relative magnitude of the vacuum (S2S) and the magnetic 
moment induced ((J,Bj_) mixings. In particular, a critical value of the ratio of the two mixing 
elements exists, 

( s ^ \ , fl + r , K1 v 
tana; c = W . (51) 



\fiB ± J c V 1 -r 

If tana; = tanu; c , the eigenvalue H^(N) of the Hamiltonian (38) will coincide with the 

flavour level energy H^(N) for all values of N (Fig. 4a): 

H 3 (N) = H^N) = -(1 + r)N + 2c 2 5. (52) 
The dependence on N of the other two eigenvalues of (38) can also be found: 



tfi,2 = § ±V4 +A (53) 
where H s is defined in eq. (10). There is a simple interpretation of condition (51) and the 
critical regime of transitions associated with it. Indeed, eq. (51) can be rewritten as 

\H e -H^\ _ (s 2 5) 2 



Ha — H e \ (/i_Bj 



Y2 



(54) 



Now, the mixing, as is well known, leads to a "repulsion" of the levels: the energy splitting 
between the levels becomes bigger than that without mixing. The larger the mixing and 
the smaller the splitting between the corresponding flavour levels, the larger the repulsion 
effect. In our case two states influence the if 3 level (which coincides with that of v^ L at 
high densities). Then eq. (54) is nothing else but a balance condition ensuring that the 
repulsion of the level from the u eL and u eR ones compensate each other for all densities. 
In other words, the shift of the H 3 level resulting from the flavour mixing is compensated 
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by that stemming from the magnetic-moment induced mixing. Due to this balance the v\ 
level does not deviate from the one and turns out to be a straight line coinciding with 



In the critical case the minimal splittings between all levels happen in the merging point. 
They are determined by eqs. (40). 

The case tana; < tanc<j c corresponds to relatively weak vacuum mixing. Now the levels 
of the neutrinos participating in the spin- flavour conversion, i.e. u e R and u^ Ll repulse each 
other stronger than the v eL and levels, and H 3 (N) will not coincide with H^N) (Fig. 
4b). The points of minimal splitting between the eigenvalues are shifted away from the 
merging point symmetrically to N m ± AN m . Moreover, the splitting itself becomes smaller 
than /. With decreasing the ratio (s2S)/(/j,B±) the shift of the minima from the merging 
point increases approaching the maximal value 



The densities N m ± AiV m correspond to the points at which the z/ e ^-level crosses the levels 
of the eigenstates of the Hamiltonian describing the {u eRl u^ L ) system when s 2 5 = (i.e. 
when the u eL -\eve\ decouples). The level splitting in these points 



is determined by the (smaller) flavour mixing. 

For tana; > tanw c the vacuum mixing dominates over the magnetic moment mixing. 
The levels of the neutrinos taking part in the spin-flavour conversion effectively attract each 
other now (Fig. 4c). Again, there are two points of minimal splitting, and they are shifted 
with respect to the merging point. When S2<V(/ i -^J-) increases the shift AN m in the points 
of minimal splitting converges to 



AN m = 



fJ.B± 



(55) 





(56) 



AN m = 



S28 



(57) 
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and the splitting in these points becomes 

AH = 2^^-iiB ± . (58) 

The shifts (57) can be deduced from the densities in the crossing points between the v eR level 
and the levels corresponding to the eigenvalues of the (f e L, v^l) Hamiltonian when fJ,B± = 
(i.e. when the v eR level decouples). The splitting is determined only by the magnetic field 
induced mixing. 

The corrections due to the state give rise to small (~ e) shifts of the levels and 
to a breaking of the symmetry of the level-crossing scheme (with respect to the reflections 
(N - N m ) - —{N - N m ), ^ - -H t ). 



6 Resonant Conversions in the Merging Region 

Let us discuss further the resonant transitions in the triple resonance region. Now three 
neutrino levels are involved simultaneously in the relevant transitions and the adiabaticity 
conditions are more complicated than in the two-neutrino case. One can introduce partial 
adiabaticity parameters i,j = 1,2,3, characterizing the "jump" probabilities between 
given pairs of levels, Hi(N) and Hj(N). The partial adiabaticity parameter can be 
defined as a ratio of the spatial width of the corresponding resonance region, Ar^, and 
the oscillation (precession) length at the resonance, l R : = nAr^/1%. The oscilla- 

tion/precession length is determined by the minimal level splitting: l 1 ^ = 2ir/AH™ n , 
AH% in = min AH^N) = min (H^N) - Hj(N)) = Hi(N™ in ) - H^N^™). The width of 
the resonance Ar^ can be estimated as the distance over which the level splitting becomes 
by a factor \[2 bigger than the minimal one: 

Arg^2|(iV) rcj ,|- 1 AiV ij -, (59) 

where AiV^ = \N — iV™ n | is determined by the condition 

lAHi^Ng™ ± AAy | = V2\AH% in \. (60) 
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This definition is equivalent to the definition considered in Section 4 for the two-neutrino 
case. 

An interesting feature of the three level system under discussion is that the is 3 level 
approaches (as function of N) the v\ and v 2 levels not symmetrically with respect to the 
point of minimal splitting N min , i.e. the way the levels approach each other when iV — > N min 
depends on whether iV < N min or N > N min . (This is equivalent to having two different 
dependences of N(t) on t for iV > N m i n and A" < N m i n (see Fig. 4)). In this case one could 
introduce two adiabaticity parameters for each pair of levels, = k^\N > N^ in ), and 
= k^(N < N% in ), as well as an average adiabaticity parameter R^ = (k^ +k^)/2. 

Consider now the resonant transitions for different values of the ratio (s 2 5)/(fj,B±). 

For critical values of the parameters (Fig. 4a) the minimal splittings for all three types 
of transitions occur at the merging point. The adiabaticity parameter determining the jump 
probability from the v\ level (which approaches u eL for sufficiently large values of N) to the 
u 3 level (~ u^l) is given by 

R^) = V2-^\(N) res \- 1 (61) 
1 — r z 

(we have used the relations AH(N) = H e - H„ for > N min , and AH(N) = H e - for 
< N m i n , and the fact that AH m i n = f, in obtaining eq. (61)). The adiabaticity parameter 
for the transition z/ 2 <-> v 3 coincides with r( 13 \ One has for the transition V\ ^ i/ 2 

K ( 12 ) = 2v / 2/ 2 |(A>) res |- 1 . (62) 

Note that all adiabaticity parameters are determined by the same splitting parameter /, 
and if r is not too close to 1, are of the same order. If the adiabaticity conditions for the 
transitions of V\ are fulfilled > 1, > 1), and v eL is produced far from the merging 

point at Ao ^> N m , it will be almost completely converted into u e R provided A" < N s f. 
If the adiabaticity is weakly broken, u e L will be transformed with large probability into 
9 e jt, with rather small probability into u^l, and with even smaller probability will remain 
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i> e L- For ~ the probabilities of the v e L v^l and v e L — > v e R transitions will be 
comparable. 

Suppose next that is produced at N ^> N rn . Then the neutrino system will be on 
the u 3 level and in the case of adiabatic change of density (i.e. > 1) will end up as 
at N = 0. However, the flavour content of the z/ 3 state (which coincides at high densities 
with u^l) changes with the variation of density: 

\u 3 >= {fjLB ± \v eL > +s 2 5\P eR > -(l + r)(N-N m )^ >} R, (63) 

where fxB± and s 2 5 obey eq. (51) and R{N) is a normalization factor. Note that the relative 
admixtures of the v e L and the P e R states are fixed, and the u^l admixture changes with N . 
In the merging point the latter is zero, i.e. the v 3 state is composed only of u e L and V e R. 

For the noncritical values of s 2 5/(fxB±) (Figs. 4b, c) the results in the case of adiabatic 
transitions are the same {y e L — > v e Ri etc.), but the types of the transitions change when 
the adiabaticity conditions do not hold. Now besides the merging point there are two other 
points of minimal splitting, at N — Ni and at iV = N 2 . The splittings and the adiabaticity 
conditions in these two points are determined by the smallest mixing (flavour or magnetic 
moment induced) rather than by /. Thus, the "jumps" are more probable in the points N 1 
and N 2 than in the merging point itself. Moreover, all conditions governing the neutrino 
transitions in the regions of Ni and N 2 are identical (up to corrections ~ e). Therefore the 
splittings, the adiabaticity parameters, and the amplitudes of the transitions are the same 
in Ni and N 2 (if, of course, \N\ has the same value in the two points). 

In the case of relatively weak vacuum mixing (s 2 S < /iB±) the u eL trajectory crosses 
both minimal splitting regions. The scheme of the transitions for v eL is similar to that in 
the critical regime: if the adiabaticity is only weakly broken, u eL will be transformed mainly 
into P e R, and with a smaller probability into (the probability of u e L to remain u e L being 
strongly suppressed). In the case of strong adiabaticity violation in the points Nip, v e L 
remains with large probability u eL and only small and equal admixtures of v e R and u^ L in 
the final state will appear. Now the two-step process u eL — > u^ L — > u eL , whose probability 
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depends on the adiabaticity in the two minimal splitting points, can be more efficient than 
the direct "jump" u e L — > v e L (from the level v\ to the level v<i) which depends on the big 
splitting 2f. 

The neutrino transitions possess interesting features also when the initial state is v^. Its 
trajectory crosses the two minimal splitting regions as well. For sufficiently weak violation of 
the adiabaticity will be transformed into u e L and u e R with rather small and comparable 
probabilities. In the case of strong adiabaticity violation will be converted mainly into 
u e R. Moreover, the stronger the violation of the adiabaticity of the transition, the smaller 
the v^l survival probability, in contrast to the behaviour of the survival probability in 
the usual two-neutrino case. 

Suppose further that the vacuum mixing is relatively strong (s 2 <5 > /i-Bj_, Fig. 4c). 
Now the conditions in the two minimal splitting points, governing the conversions of u e L 
produced at iVo N m , can be different from those in the previously considered case. 
Indeed, at N = Ni the splittings between the v x level and the other two levels turn out to 
be appreciable, so that the passing through the first minimum splitting region at N — Ni will 
not influence the is e L state considerably. Thus, neutrinos produced as v e L will cross only one 
minimal splitting region (at N = N2 < N m ) on their way from high to low density domains 
and in the case of strong adiabaticity violation in this region they will be transformed into 
u^l. The behaviour of the neutrino system in the case of V eR in the initial state can be 
formally obtained from that with v e ^ in the initial state by interchanging the configurations 
shown on Figs. 4b and 4c. 

Let us comment finally on the two mechanisms of u e L — > u e R transitions. The u e L — > v e R 
resonant conversion due to the direct (second order) u e L — u e R mixing and the two-step 
conversion involving two resonant transitions v e L — > and — > v e R coincide (are 
indistinguishable) at the merging point. One can distinguish the two different mechanisms 
of v c l — > v e R conversion if the resonance regions of the v e L — > P e R, u e L — > z^l, and — > V e R 
transitions are sufficiently well separated. 
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The direct resonant conversion takes place in one resonance region, the corresponding 
resonance condition does not depend (in the leading approximation) on the neutrino energy, 
and the energy dependence of the conversion probability follows from the energy dependence 
of the relevant adiabaticity parameter. The level splitting at the resonance point is deter- 
mined by the product of the flavour and the magnetic field induced mixing elements and 
typically is rather small. In contrast, in the case of the two-step conversion two resonant 
transitions take place in two spatial regions which are separated by macroscopic distance. 
A real propagates between the two resonant points. Both relevant resonance conditions 
depend on the neutrino energy. In each of the two resonance points the mixing is generic, 
and the splitting is generated by first order mixing. There are two adiabaticity parameters 
(and conditions) which depend on the derivative of iV in two different regions of the neutrino 
trajectory. 

Given similar conditions in the resonance regions of the three transitions v e L — > v e R, 
v eL —> and u^ L — > P eR , the two-step v eL — > u eR conversion will be more efficient 
because the adiabaticity conditions will be fulfilled better for each of the two transitions 
involved than for the direct conversion. Nevertheless, under certain conditions the direct 
conversion can be the dominant mechanism for the v e i — > v eR transitions. This can be 
the case in the situations corresponding to the energy level scheme of Fig. 2d provided 
N eS < N < N s f. This can also happen if iVo satisfies the conditions N eS < N < Nf (Fig. 
2a) and if, in addition, the spin-flavour transition u eR — > v^l is non-adiabatic (Fig. 2b). 
The direct v eL — > u eR transitions can also be efficient if v eR is produced in the initial state 
and propagates from N > N eS to iV — > 0. In all these cases the adiabaticity condition for 
the direct v e ^ — > u eR conversion can be fulfilled, in particular, due to relatively small value 
of the derivative N . 

From "microscopic" point of view there is no fundamental difference between the two 
mechanisms of v e L — > v eR conversion discussed above. In the direct case, the intermediate 
v iih{p i iR) state is virtual and both the (first-order) flavour and spin-flavour transitions which 
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result in the second-order v e ^ — v e R mixing can in principle occur at the same point. In 
the two-step conversion the intermediate u^l state is practically real, the two resonant 
conversions are separated by a macroscopic distance, and the transition probability can be 
represented as a product of the individual transition probabilities in the two resonances. If 
the u e L <-> v eR resonance point approaches the v e L <-> v^l and v eR v^l resonances, an 
enhancement of the v e i — 9 eR mixing takes place. In the merging point the direct v e i — v eR 
mixing is no longer suppressed and the direct v e L «-> u eR conversion becomes as efficient 
as any of the transitions caused by the generic first order mixings present in the effective 
Hamiltonian in eq. (3). At the same time when Nf and N s j approach N m the two resonant 
regions associated with the two transitions in the two-step v e i — > v eR conversion approach 
each other and merge; both resonances take place in the same point, thus recovering the 
case of the direct v eL — > v eR conversion. 

We have studied so far the merging of three resonances when <p = — 4c2#/(l + r ) < 0- 
The case of merging of the <-> v eL , v eL <-> v^ R) and u^ L <->• v^ R resonances (Fig. 2f), 
which can be realized for cf) = 4c25r/(l + r) > 0, can be analyzed along the same lines. 



7 Phenomenological Implications 

We would like to make a few remarks concerning the possible phenomenological implications 
of the results obtained above. Consider first the consequences for the solar neutrinos. 

1. It should be noted that magnetic field configurations in the Sun implying twisting 
fields were considered in the astrophysical literature without any reference to the solar 



neutrino problem [16|. What is really to be questioned is whether one should expect the 
magnitude of the field rotation velocity to be relevant for the evolution of the neutrino 
system. Let us show that this is really the case. As we have seen, the field rotation effects 
are important when 

~ V2G F N eff . (64) 
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Consider neutrino spin or spin-flavour precession in the convective zone (precession in the 
interior regions of the Sun is unable to account for the time variations of the solar neutrino 
flux). The magnitude of </> can be characterized by the curvature radius tq of the magnetic 
field lines, <fi ~ 1/Vo- The depth of the solar convective zone is ~ O.3i? , so that one can 
expect r to be of the order of 10% of the solar radius R Q . This gives <j) ~ 10/i? o ~ 3 x 1CT 15 
eV, whereas the r.h.s. of eq. (64) for the density near the bottom of the convective zone 
p ~ 0.16 g/cc is ~ 8 x 10 -15 eV, i.e. exactly of the same order of magnitude. It is remarkable 
that such apparently unrelated quantities as the Fermi constant, solar matter density in the 
convective zone and solar radius turn out to satisfy eq. (64). We would like to emphasize 
that one does not need any fine tuning of the parameter : for each value of <j) in the above- 
defined range and given neutrino energy E there exists a corresponding value of matter 
density in the convective zone for which the resonance condition is satisfied. 



2. As was pointed out in |T2| . [T3|| , field rotation effects can either enhance or suppress 
neutrino spin or spin-flavour precession in the Sun depending on the sign and magnitude 
of <f>. It can also change the energy dependence of the u e i survival probability and so 
distort the solar neutrino spectrum. Twisting magnetic field can give rise to some specific 
effects which could, in principle, be used for identification of the rotation. For example, the 
picture of the semiannual variations of the solar neutrino flux should crucially depend on 
the field rotation. These variations can be either enhanced or suppressed, or the minima 
of the flux can become asymmetric depending on the signs of the <fi in the northern and 
southern solar hemispheres |TJ|, 17]. The field rotation effects can also manifest themselves 
in additional time dependences of the solar neutrino flux which are not related to the time 
variation of the strength of the magnetic field (11-yr cycle). This could come about if the 
field-rotation profile in the Sun 0(r) is not constant in time. Another possible reason for 
short-time variations can be related to the rotation of the Sun and azimuthal dependence 
of the 4> profile. In this case one can expect the neutrino flux to vary in time with a period 
comparable to the solar day (~ 27 earth days). 
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3. The permutation of the resonances will lead to a change of the helicity and the flavour 
content of the neutrino final state. If the adiabaticity conditions hold, the change of the 
order of the MS W and the spin-flavour resonances will lead to the presence of a considerable 
v^l component in the neutrino flux at Earth, instead of a component. Unfortunately, 
this two possible components of the neutrino flux cannot be distinguished in the current 
and next generation experiments with solar (and supernova) neutrinos: the neutral current 
interactions of the low energy (E < 50 MeV) v^l and v^r are practically the same. 

4. The most interesting results of our study are associated with the z/ e ^ — > v e R transition. 
As we have shown, magnetic field rotation can modify drastically this transition: it can 
strongly enhance or suppress the v eL — > v eR conversion probability and, consequently, the 
v e R component of the neutrino flux from the Sun. In particular, the field rotation can induce 
direct v e i — > v e R resonant conversion, make the two-step conversion resonant, and lead to 
the phenomenon of merging of resonances. In the latter cases a strong v e i —>■ v e R conversion 
is possible. The conditions for direct isolated v e L — > v eR transitions are more restrictive 
than those for the usual RSFP. In particular, given the same magnetic field and matter 
density distributions, the adiabaticity parameter for the former is generally smaller. Thus, 
even if the RSFP conversion is efficient in the Sun, the direct v £ l — > v eR transition may not 
be efficient. However, the adiabaticity condition in the triple resonance can be fulfilled in 
the Sun even if the transition in the isolated v e L v e R resonance is non-adiabatic. This 
condition and the relevant adiabaticity parameters at the merging point depend on the 
neutrino energy. Thus, one can expect an enhancement of the v e R component of the solar 
neutrino flux in a definite energy range. A sizable flux of v eR from the Sun would be in 
contradiction with the upper limit $(^ e ) ~ (5 — 7)%$(V e ) obtained from the analyses of 
the Kamiokande II and LSD data [18|, [19|], and so some regions of values of the relevant 
parameters can be ruled out. Note, however, that these constraints are valid only if the v e i 
are produced at a density N < N s f or if the v eL — > u^ R RSFP conversion is non-adiabatic. 
This requires either small neutrino energies or large values of Am 2 ; in the first case the limit 
on the flux of solar v e obtained from the data of the Kamiokande and LSD experiments, 
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which are only sensitive to the high energy 8 B neutrinos, may not be applicable. This means 
that the upper bound on the flux of solar v e may be used to exclude some values of only 
for large enough values of Am 2 , Am 2 Z 2 x 1CT 7 eV 2 . 

The resonant v e — > v e conversion can be compatible with the existing upper bound 
on the v e flux from the Sun since the transition may not be adiabatic. It is therefore 
important to look for the v e flux in such solar neutrino experiments as SNO, Borexino and 
Super-Kamiokande which are expected to be very sensitive to z/ e 's. 

One should stress that z/ e 's can be produced by the Sun even in the non-resonant case 
when (f) — [§, [15| ; however the resonant v e — > v e conversion increases the v e flux significantly 



and also has a clear signature which allows one to distinguish this mechanism from the 
others: the flux of z/ e 's should have a peak at the energy defined by the condition sin 2a; = 1 
(s2<5 — l^B±), whereas the energy spectra of z/ e 's generated by other mechanisms (including 
non-resonant RSFP+neutrino oscillations) should be smooth. It is important to notice that 
the position and the width of the peak strongly depend on the magnetic field strength in 
the region of the resonance; therefore the energy spectrum of resonantly emitted solar z/ e 's 
should exhibit characteristic time dependence (in particular, the energy at which the v e 
spectrum achieves its maximum should vary in time). 

Obtaining experimental information on the relative fluxes and the energy spectra (and 
their time dependences) of the z/ e , v e and iy^ + u^ + u T + u T ) components of the solar neutrino 



flux will be crucial for testing the above results [20|. 

Magnetic field rotation can also be important for the transitions of the supernova neutri- 
nos. The field rotation in the supernovae can be induced by differential rotation of matter. 
As we have seen, it influences the order of the resonances which in turn will determine the 
fate of various neutrino species propagating in a collapsing star. Since neutrinos of different 
species have in general different mean energies, their transmutations can have observable 
consequences ]7[ ^TJ. This issue will be considered in detail elsewhere. 



30 



8 Conclusions 



For a system of flavour neutrinos with vacuum mixing and flavour-off-diagonal magnetic 
moments, the chirality dependent level shifts induced by the magnetic field rotation give 
rise to qualitatively new effects. The energy level scheme can be changed drastically by 
the magnetic field rotation. The location of and the distances between the resonances are 
changed. Two new types of resonances appear at nonzero densities, namely, resonances in 
the neutrino- antineutrino transitions u e L <-> u e R and <-> v^r. For a definite value of the 
magnetic field rotation velocity <fi (or of the neutrino energy) merging of three resonances in 
one point can take place. The types of resonances that merge depend on the sign of <p: for 
4> < 0, the u e L <-> i/pL, u e R <-> u^l and v e ^ <-> P 6 r resonances can occur at one point, whereas 
for > 0, the v e L <-> v^l, u e L <-> v^r and u^l <-> z/ M # resonances can merge (assuming 
c 2 5 > 0). For sufficiently large values of |0| permutation of the resonances occurs. 

The resonant v e L — > v e R transitions can proceed in two ways: i) direct v e L — * u e R 
resonant conversion induced by second order u e L — & e R mixing; ii) two-step u e L — > v e R 
conversion which proceeds via two consecutive resonant conversions u e L — > and — > 
u e R taking place in different regions separated by a macroscopic distance. Either of these 
mechanisms can dominate in the v e L — > v e R transitions, depending on the values of the 
parameters C2S, <fi, fj,B± and N . 

In the case of merging of three resonances the indicated two mechanisms of u eL — ■> u eR 
conversion also "merge" becoming indistinguishable. A complete u eL — > u eR conversion in 
one resonance region is possible in this case. 

In a medium with constant s 2 5, fiB± and having the values for which the three 
resonances merge, or the isolated v e L — > v e R resonance occurs, the oscillation-precession of 
v eL into its CP-conjugate state u e R can proceed with maximal depth. 

The effects discussed above can exist for solar neutrinos and for the neutrinos emitted 
by collapsing stars. 
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Figure Captions 

Figure 1. Neutrino energy levels versus density N for = and C2S > 0. The v e L and 
v^l levels are shown as solid lines, the z/^ R and the v eR levels as dash-dotted and dashed lines, 
respectively. Dotted lines represent the energy levels of the eigenstates of the Hamiltonian. 

Figure 2. Neutrino energy levels versus N for different values of the magnetic field 
rotation velocity 0. The notations are the same as in Fig. 1. Figures a — d correspond 
to the case of < 0, and figures e — h correspond to > 0; the consecutive figures a — d 
and e — h have been obtained by increasing the value of |0| (see the text for details). The 
v e L — > v^iR resonance is not shown in fugures a — d because of lack of space. 

Figure 3. The level-crossing scheme for isotopically neutral medium (r=l). The nota- 
tions are the same as in Fig. 1. 

Figure 4. The dependence of the neutrino energy levels on density in the case of triple 
resonance merging (0 = — 4c 2 5/(l + r) < 0) for: a) critical value of the ratio s 2 8/([iB±) 
(see eq. (51)), b) weak flavour mixing, tanc<j < tana; c , and c) strong flavour mixing, 
tancu > tancj c . The notations are the same as in Fig. 1. 
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